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ABSTRACT 

We study the thermal and kinetic Sunyaev-Zel'dovich (SZ) effect associated with 
superclusters of galaxies using the MareNostrum Universe SPH simulation. In par- 
ticular, we consider superclusters with characteristics (total mass, overdensity and 
number density of cluster members) similar to those of the Corona Borealis Super- 
cluster (CrB-SC). This paper is motivated by the detection at 33 GHz of a strong 
temperature decrement in the cosmic microwave background towards the core of this 
. supercluster (Genova-Santos et al. 2005, 2008). Multifrequency observations with VSA 

• ' and MITO suggest the existence of a thermal SZ effect component in the spectrum of 

' this cold spot, with a Comptonization parameter value of y = T.Sj^g'g x 10~^ (Battis- 

telli et al. 2006), which would account for roughly 25 per cent of the total observed 
decrement. 

From the SPH simulation, we identify nine (50 h~^ Mpc)'^ regions containing 
superclusters similar to CrB-SC, obtain the associated SZ maps and calculate the 
probability of finding such SZ signals arising from hot gas within the supercluster. 
Our results show that the Warm-Hot Intergalactic Medium (WHIM) lying in the 
' intercluster regions within the supercluster produces a thermal SZ effect much smaller 

than the observed value by MITO/ VSA. 

Neither can summing the contribution of small clusters and galaxy groups {Al < 
5 X 10^^ h~^ Mq) in the region explain the amplitude of the SZ signal. Our synthetic 
maps show peak j/-values significantly below the observations. Less than 0.3% are 
compatible at the lower end of the 1-sigma level, even when considering privileged 
orientations in which the filamentary structures are aligned along the line of sight. 
When we take into account the actual posterior distribution from the observations, 
the probability that WHIM can cause a thermal SZ signal like the one observed in 
the CrB-SC is < 1%, rising up to a 3.2% when the contribution of small clusters 
and galaxy groups is included. If the simulations provide a suitable description of the 
gas physics, then we must conclude that the thermal SZ component of the CrB spot 
most probably arises from an unknown galaxy cluster along the line of sight. On the 
other hand, the simulations also show that the kinetic SZ signal associated with the 
supercluster cannot provide an explanation for the remaining 75% of the observed 
cold spot in CrB. 

Key words: cosmic microwave background, cosmology: theory, cosmology: observa- 
tions, methods: N-body simulations 



2 /. Flores-Cacho et al. 



1 INTRODUCTION 

A detailed account of the baryons present in all known 
components of the local Universe gives a value for the 
baryon density parameter of JIb ~ (0.010 ± 0.003)/i~^ 
(see e.g. iFukueita et alJ Il998[ )). Primordial nucleosynthe- 
sis (iBurles et al. 20011) . observations of the Lyman-a for- 
est ( Ranch et al.l 19971 ) and measurements of the angu- 
lar po wer spectrum of the Cosmic Microwave Background 
(CMB, iRebolo et aLll2004 ISpergel et allbOOTl ) lead to val- 
ues two times higher for this parameter. It seems that 
there must exist a baryonic component in the local Universe 
where half of t he baryons w ould remain hidden. Theoreti- 
cal models (e.g. ICen fc Ostrikcr., 19991 ) based on gasdynam- 
ical simulations suggest that these missing baryons could 
be accounted for in a diffuse gas phase with temperatures 
10^ < T < lO^K and moderate overdensities (5 < 10 - 100), 
known as the 'warm/hot intergalactic medium' (WHIM). 
According to these simulations, the WHIM would be located 
in filaments connecting clusters of galaxies and in large-scale 
sheet-like structures. The existence of the WHIM compo- 
nent is a fairly robust conclusion against changes in the gas 
physics included in the simu lations, as feedback processes 
(see e.g. ICen fc Ostrikeill2006l ). It has been shown that even 
when considering nonequilibrium for major metal species, 
UV photoionization and galactic superwinds, the fraction of 
baryons in the form of WHIM is still of order 

There is observational evidence of this WHIM compo- 
nent, using absorption lines in the spectra of background 
sources, either in ultra violet wavelengths for the lower 
temp erature WHIM (e.g. ISavage et al. l ll998l : lNicas"tro et al.l 
20031). or in the X- ray range (first claimed detection by 



Nicastro et all (|2005l ) towards Mkn 421 although it is stil l 
under discussion; recently detected by iBuote et al.l (|2009l ) 
towards the Sculptor Wall). Apart from this possibility, the 
WHIM could also be detected via inverse Compton scatter- 
ing of CMB photons, the so-called Sunyaev-Zel'dovich (SZ) 
effect (jSunvaev fc Zeldovichll 19721 ). In this work we focus on 
this physical process and study both the thermal (tSZ) and 
kinetic (kSZ) signals associated with this gas phase. The 
SZ effect was first proposed for clusters of galaxies and has 
since been proven to be an extremely useful tool for studying 
cluster physics. When combining SZ measurements with X- 
ray and/or optical measurements it is possible to obtain an 
independent determina tion of the cosmological parameters 
l|Carlstrom et al.] |2000| ) . Recently, it has also been proven 
that SZ surveys are effec tive means of finding ne w galaxy 
clusters in blind surveys (|Staniszewski et al.ll2008l ). 

Although the SZ effect has only been observed ro- 
bustly in the direction of clusters of galaxies (for a re- 
view see iBirkinsha^J 1 19991), the aforementioned simula- 
tions bvlCen fc QstrikeJ (119 99') and other theoretical mod- 
els (|Persic et al.l Il988l . [1990 ; Boughii 19991 ) suggest that 
SZ surveys in superclusters could lead to a detection of 
the WHIM. The basic idea is that since the SZ effect is 
proportional to the line-of-sight integral of the electron 
pressure, filamentary structures of superclusters extending 
over several tens of megaparsecs could overcome the ex- 
pected low baryon overdensities and produce a measur- 
able signal wit h present-day observing facilitie s. Follow- 
ing these ideas, lAtrio-Barandela fc Miicketl (|2006l ) theoreti- 
cally computed the expected SZ anisotropics associated with 



this WHIM phase and found its dependence on gas and 
cosmological par ameters. Mo r eover , and based on numer - 
ical simulations. iDolag et al.l (1200 5l) : iHansen et al.l ll2005h: 
iHernandez-Monteagudo et al.l (|2006l ) and Roncarelli et al.l 
(|2007^ have provided coinciding predictions of the contribu- 
tion of diffuse (non-collap sed) gas to the thermal SZ power 
spectrum. In particular, iHallman et al.l (|2007l ) concluded 
that the contribution from unbound gas could be as high 
as 12% that of galaxy clusters. 

Observations of the Corona Borealis (CrB) superclus- 
ter carried o ut with the extended configuration (11 arcmin 
resolution) (IWatsori et al . 2003 ) of the Very Small Array 
(VSA) at 33 GHz (|Genova-Santos et alllioosl ) showed the 
existence of two strong and resolved negative features in a 
region where there are no known clusters of galaxies. Our 
study focuses on the strongest spot, the so-called "H-spot", 
which has a temperature decrement of AT = — 230±23 /iK. 
This region is in fact the most negative feature found in any 
of the maps produced by the VSA, which in total covered 
a sky area of 100 deg ^. Moreover, a detailed Gau ssian- 
ity study in the region (|Rubino-Martm et al.l |2006| ) finds 
a clear de viation (99.82%) at a ngular scales of multipole 
I ~ 500. iGenova-Santos et al.l (|2008l ) have recently pre- 
sented new VSA observations with an equivalent angular 
resolution of 6 arcmin that confirm the existence of the "H- 
spot" . 

The origin of the CrB spot is still unclear. Among other 
options, a possible explanation could be an SZ effect. Ob- 
servations performed with the MITO (Millimetre fc Infrared 
Testa Grigia Observatory) telescope (|De Petris et al.lll999l ) 
at 143, 214 and 272 GHz and an angular resolution of 16 ar- 
cmin showed that roughly 25% of the total signal has a ther- 
mal SZ spectral behaviour, yielding a Co mptonization pa- 
ramctcr oi y = 7.8lg;3 x 10"*^ at 68% C.L. (iBattistelh et all 
[2006. . 20071 ). Summarizing, we need to find an explanation 
for this tSZ component, but we also need to understand the 
origin of the remaining ~ 75% of the cold spot: even when 
removing the detected tSZ signal, the CrB-H decrement is 
a 3.9ct deviation with respect to Gaussianity. 

In this paper, we use numerical simulations to explore 
the possibility that the signal observed in the CrB-H spot 
originates by diffuse gas (WHIM) within the supercluster. 
Basically, we want to test whether or not the MNU sim- 
ulations are able to explain the detected tSZ signal by 
MITO /VSA in terms of the WHIM. Additionally, we test 
whether the same WHIM component can provide a relevant 
contribution to the total decrement in terms of its kinetic 
SZ effect. N-body numerical gasdynamical simulations have 
proved a very useful tool in establishing accurate theoret- 
ical predictions of the SZ signals induced ei ther by galaxy 



clust e rs or the WHIM component (see e.g. iDa Silva et al 
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2005 
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Spri n gcl et a l.ll200ll: iDo lag ct al."2005': 'Hanse n et al 
Hernandcz-Monteagudo et al. 2006; Ronc arelli et al. 
Hallman et all |2007| . |2009| ). However, these studies 



have mainly focused on the statistical properties of these 
signals (e.g. predictions for the SZ angular power spectrum, 
or the one-point probability distribution function of the full 
sky maps as a whole). In this study, we restrict our analy- 
ses to superclusters and we study the probability of finding 
spots produced by SZ signals such as those observed in CrB. 
In addition, we also provide predictions for the observability 
of those individual features in the maps produced by the the 
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next generation of instruments, such as the Planck satelhte 
jThe Planck CoUaborationI 120061 ). 

In Sj2]we briefly describe the main characteristics of the 
simulations used; in 33]a complete description of the method- 
ology applied in the study is given, and in f|4]we present a 
general statistic description of the maps. In ^ we present 
our main analysis for the particular case of the CrB-SC spot 
while in ^we provide an estimate of the observability of SZ 
signals with Planck. Finally, in !j7] we present a discussion 
and our conclusions obtained with this study. 



2 SIMULATIONS 

For the analysis reported in this paper we have used the 
MareNostrum Universe (hereafter MNU) SPH simulation, 
which is currently the largest cosmological N-body+SPH 
simulation carried out so far. It consists of 1024^ dark and 
1024'^ SPH particles in a cubic computational volume of 
500 h-i Mpc on a side.'^ This simulation was done at the 
MareNostrum supercomputer, from which it took the name, 
located at the Barcelona Supercomputer Center, Centro Na- 
cional de Supercomputacion (BSC-CNS).'^ Initial conditions 
were set up according to the so-called "Concordance A-CDM 
Model" with parameters f2m = 0.3, Qa = 0.7, ilh = 0.045, 
erg = 0.9 and Ho = 70 km s~^ Mp c~^ and a slope of n = 1 for 
the initial power spect rum (see iGottlober fc Yeped (|2007l ) 
and lYepes et all (120071 ) for more details of this simulation). 

The mass of each dark matter particle is ttidm = 
8.239 X 10^ h-i Mq and the mass of each gas particle is 
rrigas = 1.454 x 10^ h~^ M©, evolved from r edshift z = 40 
using theTREEPM-hSPH code CADGET-2 (|Springel et al.l 
l200ll : ISpringeill2005^ . The spatial force resolution is set to an 
equivalent Plummer gravitational softening of 15 h~^ kpc, 
and the SPH smoothing length is set to the 40th neighbour 
to each particle. 

We focused our work on the snapshot corresponding 
to redshift z = 0, for which complete catalogues of ob- 
jects are available, as described in the next subsection. Fi- 
nally, to investigate the dependence of the results on the 
mass and spatial resolutions of the simulations, we have 
also considered a re-simulation of the same snapshot us- 
ing identical initial conditions but with 2 x 256'' parti- 
cles, which are also equally distributed between gas and 
dark matter, with mgas = 8.065 x 10^° h~^ M© and 
m-DM = 5.397 X 10^^ h^^ M©. The same halo catalogues 
were available for this lower resolution simulation. 



2.1 Halo identification algorithms 

In order to find all structures and substructures within the 
distribution of 2 billion particles and to determine their 
properties we have used two different algorithms: a friends- 
of- friends analysis and a spherical overdensity halo- finder. 
We have started with the parallel vers ion of the hierarch i- 
cal friends-of-friends (FOF) algorithm (|Klvpin et al.lll999l ). 



^ http:/ /astro. ft. uam.cs/~marenostrum 
^ http://www.bsc.es 



We use a basic linking length of 0.17 of the mean in- 
terparticle separation to extract the FOF objects at red- 
shift z — 0. With this linking length we have identified 
more than 2 million objects with more than 20 DM par- 
ticles which closely fo llow a Sheth-Tormen mass function 
(jCottlober et al.ll2006l '). More than 4000 galaxy clusters with 
masses larger than 10^* h~^ Mq have been found. If one di- 
vides this linking length by 2" (n = 1, 3) substructures and 
in particular the centres (density peaks) of the objects can 
be identified. 

Spherical halos are identified using th e AMIGA-Halo- 
Finder (AHF) (|Knollmann fc Knebd [20091 ). an MPI paral- 
Iclizcd modification of the algorithm presented in I Gill et al.l 
((200j). The halo finder locates halos (as well as subhalos) 
as peaks in an adaptively smoothed density field. The local 
potential minima are computed for each peak and within 
spherical volumes the set of particles that are gravitation- 
ally bound to the peak is determined. If the peak contains 
more than 20 particles, then its virial (or truncation) ra- 
dius is calculated. The virial radius i?vir is the point at 
which the density profile drops below the virial overden- 
sity M(< flvir)/(47r_R^ir/ 3) = AvirPmcan, wherc Ayjr de- 
pends on the cosmology (|Kitavama fc Sutol Il996l ): in our 
case, Avir = 330. However, for subhalos this prescription is 
not appropriate. Subhalos exist in the dense environment of 
their host halo, where the density exceeds the virial one. In 
this case the density profile shows a characteristic upturn. 
This is the truncation radius of the the subhalo. We have 
checked that the number of halos found in our simulations 
is basically the same when using the FOF and AHF algo- 
rithms. 

These two catalogues (FOF and AHF) have been used 
as a tool to identify the regions of interest and to separate 
the contribution to the SZ signal of galaxy clusters from that 
of the diffuse gas, as we explain in ^ 



3 METHODOLOGY 

To characterize the statistical properties of the SZ signals as- 
sociated with superclusters of galaxies similar to the CrB, we 
proceed the following way. We first identify specific regions 
(superclusters) within the whole MNU simulation, which are 
selected to have similar characteristics (in terms of total 
mass, overdensity and number density of cluster objects) to 
those of CrB. All our analyses are focused on these partic- 
ular regions. Once these "mock CrB regions" are identified, 
we produce the associated SZ maps (both thermal and ki- 
netic components), provided that these regions are placed 
at the distance of CrB, using an average redshift f or the su- 
percluster oi z = 0.07 (|Genova-Santos et al.|[2005l ). and the 
cosmological parameters of the simulation. 

In order to isolate the contribution of the different phys- 
ical components to these maps, we produce five sets of maps, 
considering: (a) all gas particles; (b) all gas particles except 
those which are associated with clusters in the region; (c) 
only WHIM particles, defined as those gas particles hav- 
ing temperatures within the range 10^ to 10^ K; (d) only 
gas particles that belong to the clusters; and (e) only gas 
particles that belong to groups of galaxies. The last three 
sets of maps complement each other in the sense that their 
sum is equal to the map for case (a). Finally, and in order 
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to explore the dependence of the amphtude of the detected 
signals with the relative orientation of the supercluster, we 
have also produced sets of SZ maps using random orienta- 
tions of the simulation box with respect to the observer for 
all the sets of maps described above. 

All these steps are described in the following subsec- 
tions. 



3.1 Identification of superclusters similar to CrB 

The Corona Borealis sup ercluster (CrB-SC) has been mor- 
phologically described bv lSmall et al. I (|l997h as a flattened 
pancake of 100 h~^ Mpc a side with an esti mated depth of 
40 h~ ^ Mpc. According to the catalogue bv lEinasto et al.] 
(|200ll ). it includes eight clusters with total masses (includ- 
ing gas and dark matter) ranging from 1.5 x lO^^Mo to 
8.9 X IO^^^Mq, yie lding a total mass for the supercluster of 
(3-8) X IO^'^Mq (|Genova-Santos et al.|[2005l ') . However, the 
core of the supercluster, where the Abell clusters lie, extends 
over an area of only 20 h~^ Mpc a side. Hence, we decided 
to study cubic subvolumes of 50 h~^ Mpc a side, which were 
identified as similar to CrB-SC by using the two independent 
criteria described below: 

• Criterion 1: Over-densities. Assuming that the CrB 
supercluster is spherically symmetrical and all its mass is 
concentrated in the core, its average density would be of 
the order of 1.67 x 10^^ h"^ M0( h'^ Mpc)-^ The av- 
erage density of the simulation is of the order of 8.33 x 
10^" h"^ M0( h"^ Mpc)"^, calculated by dividing the to- 
tal mass (DM -I- gas) between the volume of the simulation. 
Thus, it is around 20 times smaller than the CrB average 
density. Given this fact, an overdensity criterion for the gas 
particles in the simulation can be used in order to locate 
CrB-like regions. 

We considered as candidate regions the subvolumes of 
50 h"'^ Mpc a side that contained the largest number of 
overdense gas particles, i.e. gas particles which show a den- 
sity larger than or equal to 20 times the average density of 
the simulation. These regions are guaranteed to have a sim- 
ilar overdensity to that of CrB supercluster and enough gas 
particles within to allow further filtering when building the 
different sets of maps. 

• Criterion 2: Resemblance to CrB. The number of clus- 
ters and the total mass within a given subvolume were the 
two characteristics considered in order to apply this crite- 
rion. A complete catalogue of the objects present in the sim- 
ulation and their masses was required. 

A subvolume of 50 h~^ Mpc a side was considered as 
a candidate for the study when it satisfied 1) baryonic 
mass of the order of the mass of the CrB supercluster, i.e. 
Afgas ^ 10^^ h~^ M0; and 2) the number of clusters within 
is greater or equal to 6, defining a cluster as an object with 
MTotai = Mgas -l- Mdm 5 X lo" h~^ Mq. Therefore, this 
criterion established minimum thresholds for the total mass 
of the candidate regions, the number of clusters within and 
the mass of these clusters. Note that, in order to be con- 
servative, these three thresholds lie below those of the CrB 
supercluster. 

Nine sub-volumes were finally selected and are listed 
in Table [1] All regions satisfy the two aforementioned cri- 
teria, and are centred around the largest concentration of 



Table 1. Sub- volumes considered in this study. Each corresponds 
to a cubic box of 50 Mpc a side. Columns 2 and 3 indicate 

the gas and total (dark+baryonic) mass within each subvolume. 
The last column indicates the number of clusters identified in the 
box using the AHF algorithm (see text for details). 



Subvol. 






.^^Total 


No. clusters 






Mq] 


[ Mq ] 


Al" ^ 5 X lO^-'' Mq 


001 


2.29 X 


10" 


1.53 X 10^'' 


14 


002 


2.30 X 


IQls 


1.54 X 10^'^ 


13 


003 


2.01 X 


IQls 


1.34 X 10^'^ 


11 


004 


2.94 X 


10" 


1.96 X 10^" 


27 


005 


1.90 X 


10" 


1.27 X 10^" 


12 


006 


2.45 X 


10" 


1.63 X lO^'^ 


14 


007 


2.85 X 


10" 


1.91 X 10^^ 


19 


008 


2.74 X 


10" 


1.83 X 10^" 


23 


009 


2.85 X 


10" 


1.90 X 10^'' 


21 



over-dense particles within each of them. None of the nine 
sub-volumes overlaps another. The last column in the table 
shows the total number of galaxy clusters (defined as halos 
with M ^ 5 X 10^^ h~^ Mq), which were taken from the 
AHF catalogue. 



3.2 Building SZ maps 

For each gas particle i, the simulation provides the position 
(ri), velocity (vi), temperature (Ti), density (pi) and SPH 
smoothing length (hi). Using the physical parameters for all 
the particles, we c an build the SZ maps following the pro- 
cedure described in lDa Silva et al] (|2000l ). Each gas particle 
has an associated mass profile given by mgas W{r — ri, hi), 
where r is the coordinate vector of the position in which we 
are calculating the mass profile, ri is the coordinate vector 
of the position of the i-th gas particle, hi is the smoothing 
length for that particle and W is the normalized spherically 
symmetric kerne l adopted in the sim ulations, which in our 
case IS given by (jSpringel et al.ll200lD : 



W{x,h,) = —3 




< X ^ 0.5 
0.5 < X < 1 
X > 1 



(1) 



where 2: = |r — ri\/hi. Thus, each particle occupies a sphere 
with a radius equal to its SPH smoothing length. 

The thermal component of the SZ effect (tSZ) in a cer- 
tain direction is given by the Comptonization y parameter: 



y - 



-axdl 



(2) 



where the integral is taken along the line of sight (los), or is 
the Thomson cross-section, Ue is the electron number den- 
sity, Te is the electronic temperature, fcs is the Boltzmann 
constant and nieC^ is the electron rest mass energy. The 
numerical evaluation of the integral in Eq.[2]is done by dis- 
cretizing the los integral as 



y '■ 



Aa:^^r,VK(|r„ -ri|,/iO (3) 



Mo nip 



where Aa; is the discrete cell size along the los. In this equa- 
tion, the a index runs over all the cells along the line of sight, 
the i sum runs over all the particles that contribute to the 
column of cells projected onto a pixel, is the electron 
mean molecular weight and nip is the proton mass. 
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To build the tSZ maps, we make a further approxima- 
tion, which consists in assuming a plane-parallel projection 
of the simulation boxes. At the distance of CrB, the error 
of using this approximation instead of doing the exact los 
integration is very small, as discussed below. Thus, for build- 
ing the maps we create a 3D cubic grid with a cell size of 
Aa; = 80h~^kpc a side and we use equation |3] to obtain the 
Comptonization parameter for each cell. 

A similar method is applied for the kinetic component 
of the SZ effect, given by the fo-parameter 

b = - j n^^ardl (4) 

where Vp = v ■ h is the peculiar velocity of the cluster pro- 
jected along the line of sight (n), and c is the speed of light. 
This equation is accordingly rewritten as 

^^AxVV«p,.W^(|r„-ri|,/i,) (5) 

This process was carried out for each subvolume, yielding y 
and h maps with a pixel size of Ax = SOh^^kpc. This physi- 
cal length is equivalent to an angular resolution of 1.4 arcmin 
if we place the boxes at the CrB supercluster distance (us- 
ing the redshift z = 0.07 and the cosmological parameters 
of the simulation). These maps are then degraded in resolu- 
tion, by a convolution with the appropriate Gaussian beam, 
in order to make predictions at the relevant angular resolu- 
tions. In particular, for the case of MITO/VSA, the maps 
are degraded to an angular resolution of 16'. 

Examples of SZ maps for one of the subvolumes with 
1.4' resolution are shown in Figure [TJ where it can be seen 
that the peak values match those of typical clusters of galax- 
ies, i.e. j/max ~ 10^"* and |6|max ^ 10^^. When comparing 
thes e numbers with other results from numerical simulations 
(e.g. iDa Silva et"ar.. .2000.: Springcl ct al.,.2001. '), it is impor- 
tant to take into account the effect of pix el scale, which in 
our ca se is significantly larger. Following iRoncarelli et al.l 
Figure [2] presents the distribution of the pixel val- 
ues for the whole set of nine maps, in the plane of Doppler 
|6| parameter versus y parameter. Note that in our case, we 
are integrating the SZ flux only from a 50 h^'^ Mpc region 
because we want to isolate the effect from the superclusters, 
and thus the peak of this distribution is shifted towards low 
values of y parameter. 

Figure|3]shows the effect of the beam convolution on the 
original map (Fig. [!}, for the case of VSA/MITO resolution 
(16'). The average y parameter in the maps with clusters and 
a resolution of 1.4' ranges from 0.73 x lO"'^ to 1.3 x 10"^, the 
mean value being < y >= 9.6 x 10~*. Note that this result 
is roughly an order of magnitude l ower than the averag e 
value of (y) = 2.6 x 10"'^ found by ISpringel et al.1 l|200ll ). 
The reason is again that in our case, we are only integrating 
a total depth of 50 h~^ Mpc, while in their computations 
they consider all the integrated contribution up to jz = 6. 

The approach described above makes use of two ap- 
proximations: a plane-parallel projection and approximate 
integral. To explore their validity, we have analytically ob- 
tained the profile for one single particle and we compared it 
to the one obtained using our map-mapping technique. This 
comparison yields a point-by-point difference in the profile 
of less than 5% when the particle is located in the centre of 
the map and of less than 10% when it is placed at a distance 



y Parameter — Resolutior '.^ 




Figure 1. Example of Sunyaev-Zel'dovich effect maps obtained 
from the simulations. The angular resolution of the maps is 1.4'. 
Top: thermal SZ effeet, given by the y-parameter map in loga- 
rithmic scale. Bottom: kinetic SZ effect, given by the logarithm 
of the absolute value of the b parameter. 
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Figure 2. Two-dimensional distribution of the pixel values for 
the set of nine simulated maps in supercluster regions. The two 
axes represent the y Compton parameter and the (modulus of 
the) b parameter. Contour levels enclose the regions with 50, 70, 
90 and 99.9 per cent of the total amount of pixels, respectively. 
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y Parameter - MITO/VSA RescLition (16') 
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Figure 3. Comptonization parameter (on a logarithmic scale) of 
tile same map shown in the top panel of figure [T] convolved with 
a Gaussian beam to match MITO/VSA resolution (16 arcmin). 

of ^ 4° along the diagonal of the subvolume. Both the in- 
tegrated value and the maximum value show even smaller 
differences (^ 3%) when compared to the theoretical ap- 
proach. We find that, in that case, the error in the central 
angular coordinate of a particle is smaller that the size of 
one pixel in our integration scheme. 

3.3 Subtraction of galaxy clusters contribution to 
the maps 

In order to characterize the SZ signal which is not associated 
with galaxy clusters in the subvolumes, the contribution of 
those particles which belong to clusters is eliminated from 
the los integration to produce a set of SZ maps that will 
contain only the contribution of WHIM and galaxy groups. 
In practice, two assumptions are needed to perform these 
computations. First, we have to adopt a definition of which 
halos are going to be considered as "galaxy clusters", and 
thus they will be removed from the maps. In this study, 
we shall define clusters as those objects more massive than 
5 X 10^'^ h^^ M© (in total mass). This cluster mass threshold 
does not significantly alter our results. We note that those 
halos can be easily identified in any of the two catalogues 
described in i]2.1l 

Secondly, we need to specify to what extent we are re- 
moving particles from a given halo. A reasonable assumption 
is to remove those particles lying within one virial radius 
(rvir) from the cluster. Estimates for the virial radius of the 
clusters are provided in the AHF catalogue, so we shall adopt 
those values for this purpose. Using the Tvir values adopting 
the mass limit of galaxies of M ^ 5 X 10^^ h-i Mq, clusters 
were extracted by using a spherical top-hat function; i.e. ne- 
glecting the contribution of every particle that lies within 
the virial radius of any given cluster (in the 3D grid). 

3.4 Rotations 

Gas in the WHIM phase is located mainly along filaments 
and sheet-like structures. Because the SZ effect probes the 
integral contribution of the electron pressure along the line 
of sight, we might expect certain orientations in which the 



gas in filaments is aligned along the los of the observer to 
yield large values for the y parameter, even comparable to 
that of a cluster of galaxies. In order to quantify this pos- 
sibility, after subtracting the clusters, we randomly rotated 
each subvolume 300 times, using rotations which were uni- 
formly distributed in the sphere and were defined by random 
triplets of the Euler angles. Each of the three Euler angles 
were pseudo-randomly chosen from a uniform distribution 
between 0° and 180°. Using these 301 maps (300 rotations 
plus the non-rotated original map) for each subvolume and 
each resolution, we were able to perform a statistical study 
of the orientation effects of the detected signals. 



4 STATISTICAL ANALYSIS OF THE SZ MAPS 

Here we present an overview of the overall statistical prop- 
erties of the simulated SZ maps, and how those properties 
depend on the different gas phases and the orientation of 
the supercluster. 

4.1 Analysis of the average values in the maps 

Table[2]shows the values of the ratio of total fiux (Y) in the 
tSZ maps, between the total maps and the different sepa- 
rated maps described in 321 for each of the nine selected sub- 
volumes. These numbers indicate the relative contributions 
of the different gas phases to the total (integrated) y-flux in 
supercluster regions like CrB. Roughly 73 % of the total flux 
in the maps comes from clusters with M ^ 5 x lO'^'^ h~^ M0 , 
and the other ~ 30 % comes from galaxy groups and the 
WHIM. Of this 30 per cent signal, roughly 60% is coming 
from groups, and 40% is coming from the WHIM. 

It is interesting to compare these values, which are as- 
sociated with superclusters, wit h those found in an av erage 
field. For example, in Fig. 12 in iHallman et"all (|2007l 'l it is 
shown that for an average field, the relative contribution is 
roughly two-thirds and one-third. Thus, in the case of su- 
perclusters, we find that the relative contribution of galaxy 
clusters to the average fiux is slightly larger than in an av- 
erage field, as one would expect for an overdense region. 

Finally, we also show in Table|3]the same ratios, but for 
the total SZ flux, i.e. including also the kSZ component. For 
this calculation, the kSZ component is added as computed 
in the Rayleigh-Jeans regime (i.e. ATrj = {—2y -f &)ro). 
In this case, the total flux is thus proportional to / ATrfSl. 
We note that in this case, the relative contribution of the 
WHIM component increases because at low-temperatures 
the kSZ contribution is comparable in amplitude (or even 
dominates) over the tSZ contribution. We find that the total 
SZ flux in supercluster regions is roughly explained by two- 
thirds contribution of galaxy clusters, ~ 15% galaxy groups 
and ~ 15% WHIM. 

4.2 One-point probability distribution functions 

A general statistical description of the maps can be 
given in terms of a fluctu ation analysis (see e.g. 
iRubifio-Martm fc SunvaevI |2003| 'I. which is based on the 
study of the one-point probability distribution function (1- 
PDF). This kind of analyses allows us to characterize the 
statistical properties of the sources below the confusion level 
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Table 2. Ratios of total flux {Y) in the j/-maps of the different 
gas phases, using MITO resolution. The cluster phase is defined 
by objects with M ^ 5 X lO'^^ h~^ Mq. Groups are defined as 
halos below that mass limit. The WHIM phase is defined in terms 
of gas temperatures, as 10^ < T < 10^ K. The last column shows 
the co-added contribution of WHIM and groups of galaxies. 



WITH CLUSTERS 



Subvol. 


^total 


^total 


^total 


^total 




^groups 


^WHIM 


^non — clusters 


001 


1.33 


5.14 


11.92 


3.59 


002 


1.35 


6.14 


9.37 


3.71 


003 


1.21 


9.41 


13.05 


5.47 


004 


1.41 


6.82 


6.68 


3.38 


005 


1.59 


10.06 


4.82 


3.26 


006 


1.37 


4.56 


15.96 


3.55 


007 


1.26 


6.84 


11.72 


4.32 


008 


1.43 


4.02 


13.00 


3.07 


009 


1.41 


5.79 


7.30 


3.23 


Average 


1.37 


6.53 


10.43 


3.73 



Table 3. Same as in Tabic E] but now for the ratio of total SZ 
flux in the maps, i.e. also including the kSZ component. The total 
SZ flux (S) is proportional to J ATdC, where the total tempera- 
ture decrement is computed as in the Rayleigh-Jeans regime (i.e. 
ATrj = i~2y + b)To). 



Subvol. 


'5'total 


'S'total 


■S'total 


^total 


^clusters 


5groups 


5WHIM 


■Snoii — clusters 


001 


1.41 


5.39 


8.04 


3.23 


002 


1.41 


6.65 


6.66 


3.33 


003 


1.24 


10.33 


11.14 


5.36 


004 


1.52 


8.15 


5.11 


3.14 


005 


1.83 


11.45 


3.89 


2.91 


006 


1.41 


4.77 


10.42 


3.27 


007 


1.32 


7.95 


7.79 


3.93 


008 


1.58 


4.33 


7.34 


2.72 


009 


1.85 


7.09 


3.85 


2.49 


Average 


1.51 


7.35 


7.14 


3.38 



in the maps. We now study separately the 1-PDF of the 
thermal and kinetic SZ components. 
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Figure 4. Averaged probability distribution for tSZ {y param- 
eter). These distributions are generated averaging over the nine 
subvolumes with clusters (top panel) and without clusters (bot- 
tom panel). The different colours correspond to the different res- 
olutions used, as described in the legend. We can see that after 
the cluster subtraction, the upper tail of the distribution is sig- 
nificantly reduced. 



4.2.1 tSZ 

Fig. H] shows the 1-PDF for the tSZ component, obtained 
from the j/-maps using different angular resolutions. Com- 
paring the two panels, it is clear that the upper tail of the 
distribution is significantly reduced after cluster subtraction. 
However, the tail does not disappear completely, yielding a 
spot with y parameter values within the la region of CrB 
observations. 

4.2.1.1 Rotations of the tSZ maps. We have carried 
out a statistical study to characterize the dependence of the 
tSZ signal on the different orientations of the observer with 
respect to the supercluster. Our aim is to quantify how ran- 
dom alignments of gas filaments along the line of sight affect 
the 1-PDF y-distributions. The top panel in Fig.[5]shows the 
nine 1-PDF curves after averaging over the 300 rotations in 
each case, while the bottom panel shows the 300 curves for a 
single subvolume. Note that when studying the upper tail of 
the 1-PDF for a single subvolume (bottom panel), the dis- 
persion introduced by the rotations is very small, suggesting 
that the (few) objects which are producing those signals are 
probably very close to spherical. In other words, this high-j/ 



tail is dominated by the contribution of clusters below the 
subtraction limit, or by galaxy groups. A very small number 
of possible orientations have yielded maximum values of the 
Comptonization parameter compatible at the 1 sigma level 
with the j/-signal measured in the CrB-SC. We come back 
to this point in fJS] 

4.2.1.2 Relative contribution of the WHIM to the 
total tSZ signal. We have also performed a detailed 
study of the separate contributions to the total tSZ effect 
from the different physical gas components: galaxy clusters, 
galaxy groups and WHIM. The different 1-PDF's associated 
with each case are shown in Fig. [S] now only for the angular 
resolution of MITO/VSA (16'). As expected from the den- 
sities and temperatures of the different gas phases, galaxy 
clusters are responsible of the high-y tail of the 1-PDF dis- 
tribution, while groups and the WHIM have relevant contri- 
butions at the lowest y-values. 

When comparing with the j/-value measured in CrB- 
SC, these simulations show that only compact objects (most 
probably galaxy clusters, although we cannot completely 
discard groups of galaxies) are the only component that 
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WITHOUT CLUSTERS - MITO/VSA resolution (16 arcmin) 1-PDF Distribution - k'lTO/VSA resolution 
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Figure 5. Top: Averaged (over the 301 rotations) one-point 
probability distribution functions of the y-maps for the nine in- 
dependent subvolumes after cluster subtraction. The red line cor- 
responds to the average of the nine curves. Bottom: 1-PDFs for 
the 301 maps obtained by rotations of a single subvolumc. Again, 
the red line corresponds to the average of all black lines. 

might build up a sufficiently large tSZ signal. Based on the 
probablility distributions obtained in our simulations, the 
WHIM is not a likely cause of the y-signal detected in the 
CrB-SC, as its probability distribution reaches the zero level 
at y values much smaller than those of CrB. We will inves- 
tigate this issue further in §5. 

4.2.2 kSZ 

Figure [7] shows the 1-PDF for the h maps obtained from 
the nine sub- volumes, at the same angular resolutions as in 
Figure [J] These distributions show that only clusters con- 
tribute to the high kSZ effect tails, as one would in principle 
expect given their typical bulk velocities. Considering the 
rotations of the kSZ maps does not increases the maximum 
kSZ signals which are detected in the maps. 

Finally, we have also studied the relative contribution 
of the kSZ component to the total SZ signal. As we can see 
in Fig. [S] the probability distribution of the temperature 
decrement seems almost unaffected when considering both 
tSZ and kSZ for all components except for the WHIM, in 
which case the probability distribution is slightly broadened 
with respect to the distribution considering tSZ only. This 
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Figure 6. Averaged y parameter probability distribution for the 
nine subvolumes. The different colours correspond to the different 
components found in the maps: all the objects, only clusters, only 
groups and only WHIM. The resolution of the maps used to com- 
pute these distributions is 16', the MITO/VSA resolution. The 
shaded area corresponds to the Icr interval of the observations in 
the CrB-SC and the horizontal dashed lines show the number of 
pixels that matches a given probability, showing that the WHIM 
fails to build up enough tSZ signal to be responsible for the spot 
observed in CrB-SC. 



WITH CLUSTERS 




Figure 7. Same as Fig. |3] but for the b parameter maps (kSZ 
component). 
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1-PDF Distribution - MITO/VSA resolution 
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Figure 8. Probability distribution of the temperature decre- 
ments averaged over the nine different subvolumes and using 
MITO/VSA resolution 16'. Each colour shows a different com- 
ponent of the map, as given in the legend; solid lines represent 
only the contribution to the temperature decrement of the tSZ ef- 
fect while dotted lines represent the total temperature decrement 
(considering both tSZ and kSZ). The shaded area corresponds to 
the 1(T interval of the observations carried out in the CrB-SC. The 
horizontal dashed lines indicate the equivalence between 1-PDF 
probability and the number of pixels in the maps. 

matches the expected result, as kSZ is only relevant to the 
total temperature decrement in regions with low tempera- 
tures, so neither clusters nor groups were expected to have 
a significant kSZ effect. 

In any case, when considering the overall contribution 
of the different phases, it seems that the inclusion of a 
kSZ component cannot help in increasing the total SZ sig- 
nal to the —230 /iK observed with the VSA in the CrB- 
H spot. As a reference. Fig. [8] shows as a shaded region 
the one sigma confidence interval measured by MITO/VSA 
(ATcf^ = — 2j/crB), and which accounts for 25% of the total 
observed decrement. 



5 THE CASE OF CORONA BOREALIS 

In the previous section, all the statistical analyses were fo- 
cused on the description of global properties of the SZ signals 
associated with supercluster regions. We now restrict our- 
selves to the particular case of the detection of a significant 
negative feature in the Corona Borealis supercluster (CrB- 
SC) with the combined MITO/VSA experiments. The goal 
is twofold: on the one hand, we want to probe if the MNU 
simulations are able to explain the detected tSZ signal by 
MITO/VSA in terms of gas associated with the superclus- 
ter. On the other hand, we want to explore if the kSZ signal 
associated with the supercluster can be a plausible explana- 
tion for the rest of the observed decrement in the VSA map 
(a signal of the order of —230 /iK). As pointed out above, 
even when we remove the observed tSZ component from the 
total VSA decrement, the remaining signal (roughly 75% 
of the total) is still difficult to reconcile with a primordial 
Gaussian CMB spot. 

Given that the details of the sky coverage may introduce 
biases in the results, throughout this section we restrict our 



analyses to a sky coverage equivalent to the original VSA 
sur vey in which the cold spot w as discovered, i.e. ~ 24 deg^ 
(see|G enova-Santos et aUbOOSl . for details). In practice, this 
means that all our analyses in each of the synthetic maps 
are restricted to a circle with an equivalent area of 24 deg'^, 
centred on the central pixel of each map. 

As shown in the previous section, the overall analysis in 
all the maps showed that the WHIM cannot provide a high 
enough tSZ signal to be the sole cause of the j/-signal ob- 
served in the CrB-SC. When we include the contribution of 
groups of galaxies and small galaxy clusters (below the mass 
threshold for extraction), we cannot rule out completely the 
possibility that we might indeed be facing the cause of the 
observed j/-signal. On the other hand, when considering the 
kSZ contribution, the total decrement that can be produced 
in this case is still far from the observed cold spot in the 
VSA map (-230 ^iK). 

In order to quantify how likely groups and small clus- 
ters could yield tSZ signals in the confidence interval of the 
observations, and how large the contribution of the kSZ com- 
ponent to the total signal is, we have analysed the maximum 
amplitudes obtained in the simulated maps. The top panel 
of Table |4] shows these results for the case of the full maps 
(no cluster subtraction applied), while the bottom panel of 
the same table presents the case of the maps after the clus- 
ter subtraction (hence including groups and small clusters of 
galaxies). The three columns show the maximum values for 
the y and |6| parameters, and the largest temperature decre- 
ment in the Ray leigh- Jeans regime, computed as ATrj = 
(— 2y -f b)To. The values in Tableware in good agreement 
with the typical peak amplitudes expected for the galaxy 
clusters in CrB, ba sed on their X-ray luminosi ty (see the 
estimated values in iGenova- Santos et all (120051 ) . based on 
the scaling relati on between X-ray luminosity and tSZ peak 
tempe rature by iHernandez-Monteagudo fc Rubino-MartfrJ 

diooi))- 

The values from the bottom panel of Table |4] are sig- 
nificantly smaller than those from the original maps, and 
we find that none of the subvolumes has maximum values 
within the zblcr range of the y- value observed in the CrB 
supercluster. If we consider a wider range, we find four sub- 
volumes within the ±1.3(t region of the y-signal, and all of 
them lie within the ±1.36(t interval around the Comptoniza- 
tion parameter found in the CrB-SC observations. 

That table also incorporates the maximum amplitudes 
(in absolute values) which are detected in the kSZ maps. 
While in the maps with clusters y is about one order of 
magnitude larger than h, in the maps without clusters, both 
components of the SZE are comparable in amplitude, thus 
producing a larger temperature decrement in those cases 
in which the kSZ component contributes with a negative 
sign (i.e. a clump moving away from the observer). In any 
case, the addition of the kSZ cannot provide a significant 
contribution to the total temperature decrement observed 
by the VSA (-230 ^K). 

5.1 Rotations 

Here we also investigate how random alignments of mat- 
ter along the line of sight might affect the maximum de- 
tectable SZ signal. The max;imum values of tSZ, kSZ and 
total temperature decrement in the Rayleigh-Jeans regime. 
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Table 4. Extreme values of SZ effect in the maps using 
MITO/VSA resolution. The top panel shows the results for the 
maps with clusters while the bottom panel corresponds to the 
cluster subtracted maps. Temperature values in the last column 
are obtained for the Rayleigh-Jeans region, in which wc have 
ATrj = {-2y + b)To. 



Maps with clusters 



Subvol. 


# 


max(H) 


max( 1 b| ) 


ATrj 






[10-8] 


[10-8] 


l^K] 


001 




67.4 


2.90 


-372.60 


002 




44.8 


2.40 


-244.91 


003 




82.0 


7.64 


-426.56 


004 




40.7 


5.77 


-207.81 


005 




16.4 


2.09 


-86.86 


006 




77.2 


3.13 


-418.34 


007 




87.7 


6.48 


-464.30 


008 




71.3 


2.94 


-394.63 


009 




30.9 


6.62 


-182.62 






Maps with clusters subtracted 


Subvol. 


# 


max(a) 


max{ 1 b| ) 








[10-0] 


[10-8] 




001 




1.05 


1.31 


-7.03 


002 




0.631 


1.34 


-5.64 


003 




0.613 


0.996 


-4.18 


004 




0.682 


2.24 


-2.30 


005 




0.886 


1.21 


-8.11 


006 




2.17 


1.35 


-11.60 


oor 




0.618 


0.497 


-2.50 


008 




2.01 


0.489 


-11.44 


009 




0.691 


0.763 


-4.37 



considering all 301 maps for each subvolume are summarized 
in the top section of Table [5l We also include the ranges 
of variation of these quantities (i.e. the difference between 
the highest and the lowest maximum signal found in any 
of the rotations), and the number of rotations in which the 
maps show maximum values of the Comptonization param- 
eter compatible at the 1 sigma and 1.3 sigma levels with the 
observed j/-signal in the CrB-SC. Note that when looking at 
the total temperature decrement we never obtain a value as 
high as the —230 /xK observed for the full decrement by the 
VSA. Therefore, for the rest of this subsection we focus our 
attention on the tSZ component. 

The full set of rotated y-maps also allows us to compute 
a crude estimate of the probability distribution function for 
the expected maximum values of the y parameter at this 
resolution (MITO/VSA) and sky coverage (24 deg'^), which 
is shown in Figure |9l Although the features in this curve are 
dominated by the fact that there are only 9 fully indepen- 
dent sub- volumes in the computation"^, this figure summa- 
rizes in a single plot the full range of j/- values presented in 
Tabled 

The first conclusion is that even when selecting a "priv- 
ileged" orientation of the supercluster with respect to the 
observer, the expected tSZ signals are still far from the 
MITO/VSA measurement obtained in CrB. Less than 0.3% 
of the rotations show values compatible at 1 sigma with 
the observations. Moreover, the range of variation of the 



^ The peak centered around y ~ 2.1 X 10"^ is caused by large 
groups of galaxies (with masses just below the mass threshold 
used to define clusters) in subvolumes 6 and 8. 



without clusters 
Only WHIM 
MITO/VSA Measurements 
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Figure 9. Top: Normalized probability distributions of the maxi- 
mum values of the y parameter. Black solid line shows this distri- 
bution for the maps without clusters (considering all 2709 maps; 
301 maps per each subvolume), gray solid line shows the distri- 
bution for the maps that only include WHIM. Finally, the dotted 
lino shows the posterior distribution of the SZ observations of the 
spot obtained in CrB with MITO/VSA. All distributions have 
been normalized to unity at their peak. Bottom: A zoom of the 
above plot around the region [0, 3 X 10"®] in comptonization pa- 
rameter. 



maximum signals is very small (of the order of 50% of the 
maximum signal), suggesting that the source of the high tSZ 
signals are compact and nearly spherical objects (most likely 
groups or small clusters of galaxies). 

We again emphasize that these rotations only show the 
effect of the alignment of matter that belongs to the su- 
percluster. Hence, the ranges of variation shown in table [5] 
for each subvolume are only due to gas within 50 h~^ Mpc, 
either in groups of galaxies or in a diffuse phase. 

Although the previous results firmly suggest that a y- 
signal like the one observed is unlikely caused exclusively 
by WHIM in the supercluster, we have again subjected the 
WHIM-only subvolumes to random rotations in order to ver- 
ify this statement. A summary of the results is shown in the 
bottom part of Table [S] where we can see that the maxi- 
mum tSZ effect caused by WHIM is y ^ 5.6 x lO"'^, hence, 
one order of magnitude smaller than the observations in the 
CrB-SC. The total temperature decrement also falls far from 
the observed value, always AT > —7.3 /xK, which is around 
14 times smaller than ycrB- None of the randomly rotated 
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WHIM maps yield SZ signals within the Icr interval of the 
observations. Moreover, this contribution to the total tSZ 
signal would be masked by the contribution to the tSZ line of 
sight integral by the background of the regio n, which would 
be of the order oi < y >~ 2 x 10^^ (see e.g. ISpringel et al] 
I2OOII : iRoncarelh et ai.ll2007l l . 

Hence, the second conclusion is that according to these 
simulations, the WHIM is not a likely cause of the spot, 
even at "privileged orientations", where matter is aligned 
with the observer. 

However, these conclusions do not take into account 
that the observations yielded a very broad posterior distri- 
bution, as seen in Fig.|5] Hence, we need to use a probability 
estimator that weights the probability distribution for the 
maximum values (j/max) that we obtain in the simulation 
with this broad distribution. This probability is defined as 
follows: 

P — / ^ / -fs (^max) f^ymax^ -fobs {y)dy (6) 

where Ps stands for the probability distribution for the i/max 
values in superclusters in the simulations and Pobs is the 
observed probability distribution for the spot. For this com- 
putation, we use the curves plotted in Fig. |9] (properly nor- 
malized to get a unity integral) as an estimate for Ps, while 
Pobs is assumed to have a Gaussian shape with a sigma 
of 5.3 X 10"®, taken from the MITO/VSA observations. 
The integration of Eq. |6] gives a probability of P = 3.2% 
for the maps without clusters (which still include groups) 
and a probability of P = 0.9% for the maps that only in- 
clude WHIM. Hence, supporting our conclusion that neither 
groups of galaxies, nor WHIM are likely explanations for the 
tSZ component of the spot. 



6 OBSERVABILITY OF DIFFUSE SZ SIGNALS 
FROM SUPERCLUSTERS WITH PLANCK 

Although this was not the main motivation of this paper, 
the same set of MNU simulations can be used to predict 
the expected level of SZ signals in superclusters of galaxies 
that could be detectable with the Planck satellite. To this 
end, we have convolved all the maps to a common (nominal) 
resolution of 5 arcmin and we present here a brief discussion 
of this topic by repeating the same kind of analysis that was 
done in previous section. 

We concentrate here on "diffuse" SZ signals associ- 
ated with superclusters, where for the purposes of this pa- 
per "diffuse" means "all gas phases which contribute to 
the SZ signal but excluding all clusters with masses M 
5x 10^'^ h"^ M0" (i.e. our maps with clusters subtracted). A 
detailed study of the contribution of unresolved galaxy clus- 
ters (i.e. clusters with fluxes below the detection threshold 
for Planck, but with masses greater than 5 x 10^"^ h~^ M0) 
is beyond the scope of this paper. 

Table [6] shows the expected range of variation for the 
maximum SZ signals (both tSZ and kSZ) originating in su- 
perclusters of galaxies without considering the contribution 
of the clusters themselves. The top panel includes the contri- 
bution of groups of galaxies and diffuse gas while the bottom 
panel shows only the WHIM contribution. Taking into ac- 
count the nominal sensitivities of the Planck satellite (e.g. 



6 /^K at 143 GHz with 7.1 arcmin beam, or 13 /iK at 217 GHz 
with a 5 arcmin beam), direct detections of individual fea- 
tures due to groups-!- WHIM signals might be feasible pro- 
vided that a good correction of the rest of the foreground 
components is achieved. However, a direct detection of an 
individual feature due to WHIM alone seems more difficult, 
although in principle it would be possible at the 2 or 3 sigma 
level. 

Finally, and for illustrative purposes only, we make 
here a simple comparison with the extrapolated CrB-SC y- 
parameter value, which was measured using MITO/VSA. To 
this end, we consider two reference values, namely the peak 
of the y-parameter, and the minus one-sigma value, and we 
re-scale those two numbers by assuming that scale factor 
is the one for a point-like object. In this case, we obtain 
8.0 X 10"^ and 2.6 x 10"^ for the central and minus one- 
sigma values, respectively. We note that those values are an 
order of magnitude larger than the maximum tSZ signals 
obtained in the maps without clusters. Hence, if we consider 
the Planck resolution, the numerical simulations are not able 
to explain tSZ signals as high as the one measured in CrB- 
SC (provided that the object producing it is point-like) in 
terms of hot gas in the supercluster which is not in clusters 
of galaxies. A more detailed study would require the mod- 
elling of the spot shape, which is beyond the scope of this 
article. 



7 DISCUSSION AND CONCLUSIONS 

Based on the nine supercluster regions extracted from the 
MNU simulation, our main conclusion is that a Comptoniza- 
tion y-parameter as large as the one measured in CrB-SC 
by MITO/VSA is probably produced by galaxy clusters 
(M > 5 X 10^^ h"^ Mq). When excluding the tSZ compo- 
nent associated with galaxy clusters, the maximum y-values 
are typically of the order of j/ < 2 x 10~®. Exceptionally, we 
do obtain orientations in which the maximum value is com- 
patible at the lower end of the 1 sigma level (with a largest 
value of J/ ~ 2.84 x 10~®). However, these "privileged" ori- 
entations represent less than 0.3% of all the cases. When 
we take into account the actual posterior distribution from 
the observations as described in equation [6] the probabil- 
ity of finding CrB-like features in the maps without clusters 
is 3.2%. Moreover, those values are comparable to (and in 
some cases smaller than) the expected contrib ution to the los 
SZ integral by background clusters (see e.g. ISpringel et al.l 
I2OOII : IRoncarelh et al.ll2007^ . Hence, it seems more plausi- 
ble that the j/-signal observed in the CrB-SC is caused by a 
cluster of galaxies that does not belong to the supercluster, 
and that is located further along the line of sight, instead 
of the hypothesis that the spot is being produced by small 
clusters or galaxy groups associated to the supercluster. 

On the other hand, and if our simulations provide a 
suitable description to the gas physics, the WHIM phase, 
defined as gas with 10^ K < T < lO' K, cannot build up a 
tSZ signal compatible with the observations. If we take into 
account the posterior distribution from the observations, the 
probability of finding tSZ signals compatible with the obser- 
vations (see eq. [SJ is smaller than 1%. Hence, although we 
can not rule out a contribution of WHIM from the super- 
cluster, it is unlikely the cause of the spot in CrB. 
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Table 5. Orientation effects on the extreme values of the SZ effect, using tfie resolution of MITO/VSA in the maps. The top panel 
shows the results for cluster-subtracted maps (and hence includes the contribution of both the WHIM and small groups of galaxies). The 
bottom panel shows the results for the maps only considering the WHIM contribution. The columns show: (1) subvolume identification 
number; (2) maximum value of the y parameter considering every map of the given subvolume, in units of 10~®; (3) minimum value of 
J/max obtained after rotations, in units of 10~^; (4) maximum of the absolute value of the b parameter considering every map of the given 
subvolume, in units of 10~^; (5) minimum value of |6|max obtained after rotations, in units of 10~^ (6) largest temperature decrement 
ATjnin, expressed in fj,K in the Rayleigh-Jeans regime, considering every map of the given subvolume. N.B. Temperature decrements 
are negative quantities, hence shown here are the minimum values obtained; (7) minimum temperature decrement in the Rayleigh-Jeans 
regime obtained after rotations, expressed in fiK; (8) number of independent maps that yield maximum values of the y parameter in the 
itlsigma range of CrB. (9) number of independent maps that yield maximum values of the y parameter in the ±1.3sigma range of CrB. 



Maps with clusters subtracted (Only galaxy groups and WHIM) 



Subvol. 


max(j/n,ax) 




max{ 1 &|max 


) min{|61„ax) 


min(AT,-„i„) 


max(AT,-„i„) 




N±1.3, 




[10-8] 


[10-6] 


[10-6] 


[10-6] 










001 


1.39 


0.77 


1.80 


0.68 


-11.05 


-5.64 





193 


002 


1.30 


0.42 


2.08 


0.83 


-8.83 


-3.62 





19 


003 


1.08 


0.54 


1.48 


0.91 


-6.23 


-3.22 





7 


004 


1.04 


0.58 


2.65 


1.46 


-3.80 


-1.56 





19 


006 


1.07 


0.49 


1.60 


0.91 


-10.22 


-5.44 





8 


006 


2.62 


1.25 


1.35 


0.42 


-14.63 


-6.84 


5 


301 


007 


1.01 


0.44 


1.89 


0.44 


-5.00 


-1.76 





18 


008 


2.84 


1.57 


1.03 


0.40 


-17.22 


-8.94 


3 


301 


009 


1.05 


0.60 


2.57 


0.62 


-12.72 


-4.17 





4 










Only WHIM 










Subvol. 


max(i/niax) 


min(j/niax) 


max( 1 &|max 


) min{j6|„,ax) 


min(AT,„i„) 


max(AT,„i„) 


N±i„ 


N±1.3„ 




[10-6] 


[10-6] 


[10-6] 


[10-6] 


[mX] 








001 


0.51 


0.22 


0.91 


0.36 


-3.91 


-1.51 








002 


0.53 


0.19 


1.57 


0.75 


-6.25 


-3.11 








003 


0.50 


0.25 


1.49 


0.64 


-4.55 


-2.01 








004 


0.55 


0.21 


2.25 


1.07 


-2.39 


-0.93 








006 


0.49 


0.24 


1.10 


0.49 


-4.84 


-3.01 








006 


0.42 


0.15 


1.15 


0.46 


-4.05 


-0.82 








007 


0.36 


0.15 


1.19 


0.31 


-3.28 


-1.01 








008 


0.46 


0.21 


0.82 


0.29 


-3.49 


-1.38 








009 


0.48 


0.10 


1.73 


0.43 


-7.32 


-1.61 









Table 6. Forecast of the maximum expected values for the diffuse SZ signals in supercluster of galaxies for Planck. For these computations, 
we adopt an angular resolution of 5 arcmin. For the meaning of the different columns, see Table [5] for comparison. 





Maps 


with clusters 


subtracted (Only galaxy groups 


and WHIM) 




Subvol. 


max(j/n,ax) 


min{j;„ax) 


max( 6 max) mi 


in(|6|max) 


min(ATmi„) 


max(ATmin) 




[10-6] 


[10-6] 


[10-6] 


[10-6] 


[mA'] 


[mA] 


001 


3.07 


1.56 


7.25 


3.90 


-36.50 


-19.67 


002 


2.91 


0.94 


6.05 


2.98 


-26.48 


-13.05 


003 


4.26 


2.02 


4.62 


3.68 


-19.52 


-9.82 


004 


3.36 


1.42 


8.89 


4.69 


-13.85 


-4.18 


OOS 


3.55 


2.32 


5.78 


4.14 


-34.73 


-24.03 


006 


4.31 


1.74 


3.26 


1.41 


-23.17 


-10.86 


007 


2.35 


0.84 


5.67 


1.36 


-15.47 


-5.88 


008 


5.77 


2.46 


2.71 


1.49 


-30.84 


-14.45 


009 


2.87 


0.86 


5.57 


1.65 


-25.22 


-7.26 








Only WHIM 








Subvol. 


max(j/„ax) 


min{j;,„ax) 


max( 6 max) nri 


in(|b|max) 


min(ATmi„) 


max(ATmin) 




[10-6] 


[10-6] 


[10-6] 


[10-6] 


[mA] 


[mA] 


001 


1.94 


0.80 


3.26 


1.57 


-14.77 


-5.99 


002 


1.76 


0.86 


6.06 


3.22 


-25.99 


-12.93 


003 


1.64 


0.83 


3.91 


1.78 


-16.20 


-8.18 


004 


1.69 


0.88 


7.69 


3.54 


-8.37 


-3.85 


005 


1.55 


0.67 


2.70 


1.38 


-12.47 


-7.05 


006 


1.35 


0.49 


2.72 


1.46 


-11.14 


-2.80 


007 


1.14 


0.69 


3.79 


1.08 


-9.32 


-5.51 


008 


1.37 


0.84 


2.52 


1.41 


-10.88 


-6.74 


009 


1.31 


0.44 


5.32 


1.59 


-18.85 


-6.69 
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Including the kSZ component does not modify this con- 
clusion. However, the simulations show that the kinematic 
component of the SZE plays an important role in building 
up the total temperature decrement in the maps with clus- 
ters removed. Whereas in clusters kSZ is about one order of 
magnitude lower than tSZ, when clusters are excluded from 
the analysis, both SZ components are comparable. This role 
is also noticeable in the temperature decrement probabil- 
ity distribution, showing that in order to study the SZ effect 
due to WHIM, kSZ should not be neg lected. This idea seems 
to be supported by the work done bv lAtrio-Barandela et all 
(|2008l ). where it is shown that the kSZ effect originates in 
filaments of gas of certain characteristics could yield tem- 
perature decrements comparable to those of CrB. However, 
it should be noted that we have not been able to repro- 
duce the temperature anisotropics they predict, considering 
galaxy groups plus WHIM in our simulations. 

Our main conclusion might in principle depend on the 
cosmology. The MNU simulation has been obtained using 
a rather large normalization power spectrum to set up the 
initial conditions. Current WMAP5 results indicate smaller 
values for both og, and fim than the ones we used in the 
simulation. This will translate into less evolved structures, 
where the abundance of massive clusters in t he same vol- 
ume is considerably suppressed as was shown in lYepes et al.l 
(|2007l ) from a set of diflerent MareNostrum simulations with 
WMAPl and WMAP3 cosmological parameters. This means 
that, attending only to the cosmology, the results shown in 
this paper constitute an upper limit to the probability of 
finding features like the CrB spot for the same cosmic vol- 
ume sampled. We cannot neglect either the effect of cosmic 
variance: using a larger volume we would have been able 
to extend the sample of superclusters and thus improve the 
statistics in order to obtain a better probability function. 
However, according to the variance shown in the simulation, 
we do not expect significant changes in the results. 

The other caveat concerns the modelling of the physics 
of baryons in clusters. MNU is a pure adiabatic gasdynam- 
ical simulation where radiative processes are not taken into 
account. The inclusion of radiative cooling, UV photoion- 
ization and star formation would change the properties of 
the gas inside clusters and the corresponding SZ effect. How 
large this change is not a simple answer since we are not 
yet able to perform large scale simulations of the size of 
the MareNostrum Universe but inc luding also cooling an d 
star formation effects. According to lCen fc Ostrikeij (|2006t ). 
the inclusion of these feedback effects does not significantly 
modify the fraction of baryons expected in the WHIM. How- 
ever, feedback does slightly modify the differential mass frac- 
tion, heating the gas from the warm phase (T < KT" K) 
into the WHIM phase, while the hot phase (T > lO' K) 
remains unchanged. Given that the greatest SZ effect from 
the WHIM originates in the phas e with T ~ lO' K (e.g. 
iHernandez-Monteagudo et al.ll2006l ). our results seem fairly 
robust against feedback. Nonetheless, it should be further 
studied, especially since shocks could play a very important 
role in heating the gas in the CrB cluster and hence be the 
cause of the spot. Despite the lack of large scale simulations 
with feedback effects, one possibility for this study is to per- 
form a re-simulation of a volume-limited sample of clusters 
in which considerably more physical processes are to be in- 
cluded. This will allow us to estimate the relative change in 



the determination of the SZ effect with respect to the adia- 
batic runs and estimate the effect of shocks. This is a future 
project. 

All in all, and according to the description provided by 
the MNU gasdynamical simulations, we conclude that the 
observed thermal SZ component of the CrB H-spot most 
probably arises from one (or several) unknown galaxy clus- 
ter(s) along the line of sight, which in principle are not 
neccesarily associated with the CrB supercluster. In order 
to explore this idea further, we have started an observa- 
tional programme to perform a detailed characterization of 
the properties of the gala xy populations along the line of 
sight of the CrB spot (see I Padilla- Torres et al]|2009l '). 
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